Numerical integration of the equations of motion for rigid polyatomics: The matrix method
نویسنده
چکیده
A new scheme for numerical integration of motion for classical systems composed of rigid polyatomic molecules is proposed. The scheme is based on a matrix representation of the rotational degrees of freedom. The equations of motion are integrated within the Verlet framework in velocity form. It is shown that, contrary to previous methods, in the approach introduced the rigidity of molecules can be conserved automatically without any additional transformations. A comparison of various techniques with respect to numerical stability is made.
منابع مشابه
On the numerical integration of motion for rigid polyatomics: The modified quaternion approach
A revised version of the quaternion approach for numerical integration of the equations of motion for rigid polyatomic molecules is proposed. The modified approach is based on a formulation of the quaternion dynamics with constraints. This allows to resolve the rigidity problem rigorously using constraint forces. It is shown that the procedure for preservation of molecular rigidity can be reali...
متن کاملAn Analytical Technique for Solving Nonlinear Oscillators of the Motion of a Rigid Rod Rocking Bock and Tapered Beams
In this paper, a new analytical approach has been presented for solving strongly nonlinear oscillator problems. Iteration perturbation method leads us to high accurate solution. Two different high nonlinear examples are also presented to show the application and accuracy of the presented method. The results are compared with analytical methods and with the numerical solution using Runge-Kutta m...
متن کاملAn Efficient Numerical Algorithm For Solving Linear Differential Equations of Arbitrary Order And Coefficients
Referring to one of the recent works of the authors, presented in~cite{differentialbpf}, for numerical solution of linear differential equations, an alternative scheme is proposed in this article to considerably improve the accuracy and efficiency. For this purpose, triangular functions as a set of orthogonal functions are used. By using a special representation of the vector forms of triangula...
متن کاملA new approach for solving the first-order linear matrix differential equations
Abstract. The main contribution of the current paper is to propose a new effective numerical method for solving the first-order linear matrix differential equations. Properties of the Legendre basis operational matrix of integration together with a collocation method are applied to reduce the problem to a coupled linear matrix equations. Afterwards, an iterative algorithm is examined for solvin...
متن کاملA computational wavelet method for numerical solution of stochastic Volterra-Fredholm integral equations
A Legendre wavelet method is presented for numerical solutions of stochastic Volterra-Fredholm integral equations. The main characteristic of the proposed method is that it reduces stochastic Volterra-Fredholm integral equations into a linear system of equations. Convergence and error analysis of the Legendre wavelets basis are investigated. The efficiency and accuracy of the proposed method wa...
متن کامل